THRESHOLDS FOR SHOCK FORMATION IN TRAFFIC FLOW 
MODELS WITH ARRHENIUS LOOK-AHEAD DYNAMICS 

YONGKI LEE AND HAILIANG LIU 

Abstract. We investigate a class of nonlocal conservation laws with the nonlinear ad- 

vcction coupling both local and nonlocal mechanism, which arises in several applications 

CnI ■ such as the collective motion of cells and traffic flows. It is proved that the C 1 solution 

Vh ' regularity of this class of conservation laws will persist at least for a short time. This 

Mh. persistency may continue as long as the solution gradient remains bounded. Based on 

■^^ \ this result, we further identify sub-thresholds for finite time shock formation in traffic 

flow models with Arrhcnius look-ahead dynamics. 



m 



1. Introduction 



< 

In this work we investigate a class of nonlocal conservation laws, 

, v f d t u + d x F(u, u) = 0, t>0,xGl, 

^ ' ' \ u(0,x) = u (x), iGK, 

>• ! where u is the unknown, F is a given smooth function, and u is given by 

^: (1.2) u(t,x) = (K*u)(t,x) = I K(x-y)u(t,y)dy, 



•*^£ \ where K is assumed in PV" 1 ' 1 (R). The advection couples both local and nonlocal mech- 

anism. This class of conservation laws appears in several applications including traffic 
flows [HI ED], the collective motion of biological cells [31 EH CLT], dispersive water waves 
[23| El [21 [10] , the radiating gas motion [51 [191 [IS] an d high-frequency waves in relaxing 
medium [H[2T1[22]. 

r^ We are interested in the persistence of the C l solution regularity for (II .11) . As is known 

that the typical well-posedness result asserts that either a solution of a time-dependent 
PDE exists for all time or else there is a finite time such that some norm of the solution 
becomes unbounded as the life span is approached. The natural question is whether there 
is a critical threshold for the initial data such that the persistence of the C 1 solution 
regularity depends only on crossing such a critical threshold. This concept of critical 
threshold and associated methodology is originated and developed in a series of papers 
by Engelberg, Liu and Tadmor [H [121 C3] for a class of Euler-Poisson equations. 

In this paper we attempt to study such a critical phenomena in (11.11) . C 1 solution 
regularity is shown to persist at least for finite time. Moreover, such persistency may 
continue as long as the solution gradient remains bounded. We also identify sub-thresholds 
for finite time shock formation in some special traffic flow models, as well as fll.l[) with 
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one sided interaction kernels. These together partially confirm the critical threshold 
phenomenon in non-local conservation laws (11.11) . 

The traffic flow model that motivated this study is the one with looking ahead relaxation 
introduced by Sopasakis and Katsoulakis [20] : 

d t u + dju(l - u)e~ K * u ) = 0, t > 0, x G K, 



u(0, x) = Uq{x), X G 

where u(t, x) represents a vehicle density normalized in the interval [0, 1] and the relaxation 
kernel 



(1.4) K(r) 



7 

0, otherwise 



if _ 7 < r < o, 



is the constant interaction potential, where 7 is a positive constant proportional to the 
look-ahead distance and K Q is a positive interaction strength. We set K = 1 since in our 
study this parameter is not essential. 

An improved interaction potential for (jl.3p is introduced in [8] with 

(1.5) K(r) = (f(l + ?). -^r<0, 

(^ 0, otherwise. 

This linear potential is intended to take into account the fact that a car's speed is affected 
more by nearby vehicles than distant ones. The authors in [8] carried out some careful 
numerical study of the traffic flow model JO]), through three examples: red light traffic, 
traffic jam on a busy freeway and a numerical breakdown study. In the case of a good 
visibility (large 7), their numerical studies suggest that ( 11. 3ft with the modified potential 
( II. 5p yields solutions that seem to better correspond to reality. 

The objective of this article is therefore twofold : i) to establish local wellposedness of 
smooth solutions for (II. ip ; ii) to identify threshold conditions for the finite time shock 
formation of the traffic flow model (II. 3p subject to two different potentials (II .4p and ( II. 5p . 
respectively. The finite time shock formation of solutions in traffic flows are understood 
as congestion formation. 

We use X to denote a space X(M) for X = H 2 , W 1 ' 1 and L°° + H 2 . The main results 
are collectively stated as follows. 

Theorem 1.1. (Local existence) Suppose F G C 3 (1R, M.) and K G W 1 ' 1 . If Uq G 

L°° + H 2 , then there exists T > 0, depending on the data, such that (ll.ip admits a unique 
solution u G C 1 ([0,T) x R). Moreover, if the maximum life span T* < 00, then 

lim \\d x u(t, -)\\ L oo = 00. 

Theorem 1.2. Consider (jl.3p with constant potential (II. 4p . Suppose that u G H 2 and 
< u (x) < 1 for all x G R. If 

(1.6) sup[u' Q (x)} > - ( - + — • 3 - min { - 1, 7 • inf [u' (x)]} 
then u x must blow up at some finite time. 
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Theorem 1.3. Consider ( 11. 3p with linear potential (ll.5p . Suppose that uq G -ff 2 and 
< u (x) < 1 /or all x G M. 1/ 

(1.7) supK(x)] > -(1 + -- /6-min{ -2, 7 • inf [u(,(x)]} 

£/*en u x must blow up at some finite time. 



Regarding these results several remarks are in order. 

i) Our threshold results in Theorems 11.21 and 11.31 are valid for any < 7 < 00. When the 
look-ahead distance 7 — > 00, both threshold conditions are reduced to sup xGR [1*0(0;)] > 0. 
On the other hand, when 7 — > 00, model (jl.3p is reduced to the classical Lightwill- 
Whitham-Richards(LWR) model [H E 



dtu + d x (u(l-u)) = 0. 

This local model can be verified to have finite time shock formation if initial data has 
positive slope i* > at some point. Therefore, the threshold conditions identified are 
consistent with that of the LWR model. 

ii) In a recent work [9] D. Li and T. Li presents several finite time shock formation 
scenarios of solutions to (11.31) with (ll.4p . Their approach is to analyze the solutions along 
two characteristic lines defined by = u(t,Xi(t)) and 1 = it(t, ^(t)), with which they 
justified that if there exist two points a.\ < a 2 , such that 1*0(0:1) = and u (a 2 ) = 1, 
then u x must blow up at some finite time. Compare to their result, our shock formation 
conditions in Theorems 11.21 and ll .31 may be viewed in the perspective of critical thresholds. 

iii) The shock formation conditions in Theorem 11.21 and 11.31 are consistent with the nu- 
merical results obtained in j5]. Indeed, a numerical comparison in [8] of solutions to (11.31) 
with (ll.4p for 7 = 0.1 and 7 = 1 indicates that the solution with 7 = 0.1 remains smooth, 
while the solution with 7 = 1 seems to contain a shock discontinuity. 

iv) The threshold in (11.71) is bigger than that in (11.61) . This observation suggests that 
under certain initial configuration, the traffic flow model with constant interaction po- 
tential may develop a congestion formation, while the model with the linear interaction 
potential may not. Roughly speaking, it is understood that the drivers with the linear 
potential are 'smarter' than the drivers with the constant potential. 

v) For fixed 7 > 0, both (11.61) and (11.71) reflect some balance between sup,j. eR [it (a;)] and 
inf x£ ]R[ii (o;)] for the finite time shock formation: if the non-positive term inf xG ]R[ito(o;)] is 
relatively small, then sup xeR [i* (x)] needs to be large for the finite time shock formation. 
It indicates that not only the car density behind the traffic jam but also the car density 
ahead of the traffic jam contribute to the formation of congestion. 
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We now summarize the main arguments in our proofs to follow. For the proof of 
Theorem II. 1\ we apply the Banach fixed-point theorem to the transformation S defined 
through v = S(u), where v is solved from 

/, on / d t v + F u v x + F u u x = 0, 

( L8 ) \v(t = 0)=u . 

We show that there exists T > depending on initial data such that the mapping v = S(u) 
exists and is a contraction. In so showing, detailed estimates of non-local terms are crucial, 
and allow us to track the dependence of T on the initial data. 

For the proofs of Theorem II .2H1 .3| we trace the Lagrangian dynamics of d := u x , which 
can be obtained from the Eulerian formulation: 

(1.9) (dt + (1 - 2u)e- n d x )d = e~ u [2d 2 + 2(1 - 2u)u x d - «(1 - u){u x } 2 + u{l - u)u xx ] . 

The right hand side is quadratic in d, the a priori bound < u < 1 ensures the bounded- 
ness of both u and u x involved in the coefficients. The key in our approach is to bound the 
nonlocal term u xx in terms of M — sup x£R [u x (x, t)] and iV = mi xeR [u x (x,t)]. This way 
we are able to obtain weakly coupled differential inequalities for both M and N, which 
yield the desired sub-thresholds. 

From the proofs of Theorem 11.2111.31 we observe that the one-sided interaction property 
of kernels (11.41) and (11.51) is crucial. Hence our threshold analysis for the traffic flow 
models is applicable to the class of nonlocal conservation laws (II. ip under the following 
assumptions: 
(ifi). F G C 3 (M,M), and the kernel K{r) e W 1 ' 1 satisfying 



Kir) 



N ondecr easing , r < 0, 
0, r > 0. 



(H 2 ). F(0,-) =F(m,-) =0 and 

F uu < 0, Fun > 0, F u < for ue[Q,m]. 

The result can be stated as follows. 

Theorem 1.4. Consider (11.11) with (11.21) under assumptions (H\)-(H.i). If Uq G H 2 and 
< uq(x) < m for all i6K, then there exists a non-increasing function A(-) such that if 

sup[«d(x)] > A(inf [u' Q (x)]), 

xeR X ^ R 

then u x must blow up at some finite time. 

We should point out that it was the threshold analysis for traffic flow models that led 
us to the thresholds (11.61) . (11.71) in the first place, which in turn was then extended to the 
general class (II. ip as summarized in Theorem 11.41 

We now conclude this section by outlining the rest of the paper. In section 2, we prove 
local wellposedness for the class of nonlocal conservation laws (11.11) . In section 3, we 
investigate sub-thresholds for nonlocal traffic flow models. We finally sketch the proof of 
Theorem 11.41 in the end of this paper. 
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2. Local wellposedness and regularity 

In this section, we study the local well-posedness of (jl.ip . We consider a solution space 
as u G uq(x) + I? T , with B T := L°°([0, T]; H%), which allows u to be non-zero at far field. 
By transformation 

U = u — u , 

we find the following equation for U G B T , 

U t + d x F(U + u ,U + n ) = 0. 

This lies in the same class as fll.ip . With this in mind, from now on, we shall consider 

ueB T :=L°°([0,T];H 2 x ). 

We prove the local wellposedness result by the fixed point argument. That is, we first 
define a transformation S as v = S(u), where v is solved from the following equation 

(0 -.\ j d t v + F u v x + F n u x = 0, 

1 j \v(t = 0)=u , 

and then show this mapping has a fixed point. 

We begin by verifying the existence of v — S(u), which is carried out in a series of 
Lemmata 2.1-2.3. For simplicity, we take 

a = F u and b = —F n u x . 

We bound a and b in terms of u in the following lemma. 

Lemma 2.1. Suppose u G B T , K G W 1 ' 1 . Then 

(2.2) \\a x \\ H i < (fc(l + ||^|| L i)) 2 (l + IKIUIMI^ 
and 

(2.3) \\b\\ H 2 < fc(l + ||K|| L i) 3 (l + \\K x \\ L i){l + ||^||oo) 2 |kl|^, 

where k = k(F) is a constant depending on F. In particular, z/sup ie r 0T i ||w||/r 2 < R> then 

sup ||a x ||//i < c a R 2 and sup ||&||h 2 < Q>-R 3 , 

te[o,T] te[o,T] 

where c a = k{\ + Ci)(l + ||if||i,i) 2 , q, = k(l + Ci) 2 (l + ||if||wi,i) 4 and c 1 is an embedding 
constant. 

Proof. We begin with some key inequalities for u: using \\w * K\\ L 2 < \\K\\ l i\\w\\l^ and 
K G W 1 ' 1 we obtain 

ll^xlU 2 = \\K*u x \\ L 2 < ||A'|| l i||m x .|| L 2, 

(2.4) llw^H^ = \\K*u xx \\ L 2 < \\K\\ L i\\u xx \\ L 2, 

\\u xxx \\ L 2 = \\K x *u xx \\ L 2 < \\K x \\ L i\\u xx \\ L 2 

and 

Halloo < H^IUII^IU 1 - 

We calculate 

"X * UU">X * " UU^XI 

&xx -f uuu'U'x ' uuuV'x'U'X \ " uu^xx \ " uuu^x^x "T "uuu^ x "T " uu^xxi 
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SO that 



iiaxiU 2 < ^ikiil 2 + ^ii^iU 1 n^iU 2 

<k(l + \\K\\ L i)\\u\\h2. 



\\o>xx\\l 2 < kl ||Mx||oo||Mx||l 2 + |kx||oo||-ft'||L 1 ||Wx||L 2 + ll^xxlli 2 

(2 5) 

+ |K-||oo||A|| L i||u x || L 2 + |hx-||oo||-ft'||ii||«x|U 2 + ||A'|| L i||u xx || L 2 

< fc(i + ||m x || 00 )(i + HA'iiii) 2 !^!!^. 

These together lead to (I272|) . 
We also calculate, 

®x -f uu^xUx "uvV'x " U^XXl 

®xx ~ " VlUu^x^x " uuu"'x'^x uu" j xU'x "ra%"ii 

-f UUU"'X"'x uuu"-x uu" J x" J xx 

" uu^x^xx " uu" , x"'xx " U^XXXl 

to obtain 

IHU 2 < fc||-^||Li||^||L 2 - 

HMU 2 < ^||mx||oo||A'|| L i||m z || L 2 + A; 1 1 w^. 1 1 oo j | JiT 1 1 ^i 1 1 w^. 1 1 i2 + £;||A"|| L i||w x . x .|| L 2 

< fc( (i + ||m x || 00 )(i + ||a|| L i) 2 ) |H| ff 2. 



I^ixIU 2 < k[ \\UxWIoWKWl1 + HWillLll-^llii + IK||oo||A|| L i + \\u x \\ \\K\\ L i )IM|tf 2 



+ kn\u x \\t \\K\\ z Ll + II^II^HAlli! + 2||« a! || 00 ||A r ||ii I ||«||Ha 

+ kl lluxllooll-RTHii + ||Wx||tx>||-K"||i,i + II^-IIl 1 J \\u\\h 2 - 

These estimates give (j2.3p . D 

Lemma 2.2 (A priori estimates). Suppose u e £? T . ,4 sufficiently smooth solution v of 
(12. ip mus£ satisfy the energy estimates 

(2.6) sup \\v(-,t)\\ L 2 < (\\uo\\ L 2+T- sup ||6|| L 2 ) exp I - / ||ax||oo^ 



te[o,T] V te[o,T] J \2 Jo 



(2.7) sup \\v{-,t)\\ H 2 < ( ||n ||//2 +T- sup ||6||#2 ) exp ( (- + a) \ 
te[o,T] V te[o,T\ J \ z Jo 

where C\ is an embedding constant. 

Proof. Apply d l x to the first equation of (12.11) to obtain, 

(2.8) (d l x v) t + a-(d l x v) x = h l , 
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where h = d l x b — d l x (av x ) + a{cr x v) x . Multiplying (12. 8 j) by d l x v and integrating over R, we 
obtain, 

(2-9) ~ l{d l x vfdx= fa x { ^ + fh l -{tf x v)dx. 

z al JR JR z JR 

This with I = leads to 
d 

It 
That is 



,|li, = / 0l ^ + 2/^<| k |U|,& + 2 Wi ,|HU, 



— IklU 2 < oll a x-lloo||w||L2 + IHU 2 , 



d „ 1 

— m r2 < — 

dt" " - 2 1 
which upon integration gives (12.6J1 . Next, summing (12.91) for / = 0, 1, 2, we obtain 

(2-10) 



1 rf II ||2 

2Jt Mm 2 



\f a x -^{d l x v?dx+ f^-^dx 

JR l=Q JR l=Q 

- / a x (v 2 -v 2 x - 3vl x ) dx - a xx v x v xx dx + / (few + b x v x + b xx )v xx dx 

^ JR JR JR 



3 

< — 1 1 a^- 1 1 oo 1 1 ^ 1 1 1/2 + ||wx-||oo||a^||L2|K ; ^IU 2 + H^llff 2 \\v\\m 



Therefore, we obtain 



-jdMlfp < ( o +Cl 



which upon integration again gives (12. 7p . D 

Lemma 2.3. Suppose the initial data v(x, 0) =u 6 -f^ 2 - T/ien /or eac/i w G -B T , i/iere 
exists a unique solution v G 5 T o/ (12. ip . 

Proof. Since sup tG[or] ||o x ||hi < oo, 

— - = a, x(0) = x 
at 

admits a unique solution a; = x(x ,t) for each x G R. Along x(x ,t), (12. ip reduces to 

dv 

— = b, u(0) =M (x ). 

Hence w(x(xo, £), £) = Uo(#o) + J b(x(xo, t), r) rfr and the unique solution for (12. ip exists. 

D 

Now, we are ready to prove Theorem 11.11 

Proof of Theorem 11.11 Let R be any number satisfying R > 2||w ||#2, we define 

(2.11) B T R := iwGL°°([0,T];# 2 )|u;(x,0) = Wo, sup ||w(-, t)\\ H * < r\. 

I te[o,T) J 



YONGKI LEE AND HAILIANG LIU 



Assume that u e B^, we then have 

|K*)||oo<Co#, ||«x(t)||oo<Cii2, 0<£<T, 

where Co and Ci are the embedding constants. 

We first show that S maps B^ into B^ for some T small. From (j2.7p . it follows that 

sup |K-,£)||^< (^ + T-c 6 i? 3 )exp(V(~ + Cl ) C J? 2 ) 

(2.12) te[o,T] \2 / V l J 

<R, 

provided 

1 
" i:= 3(2 + Cl )(c a + c b )eR 2 ' 
Hence, 

S :Bl^ B T R , VT < T\. 

We next show that S is a contraction on 5^ in the L°°([0, T]; L^) norm: 

1 
te[o,T] 2 t g[o,T] 



(2.13) sup HS'(wi) - S(u 2 )\\l 2 < - ■ sup ||iti — W2IU2, Vwi,w 2 e .B^. 



Let u := wi — W2 = S'(wi) — S'(m2), then difference of ( 12. ip for t> 2 and i>i, respectively, leads 
to 

(2.14) f; t + o(mi)t; x = b, v(0, x) = 
with 

(2.15) 6 =: -{o(ui) - a(n 2 )}f 2:c + 6( w - H^)- 
Applying ( 12. 6 p we have 

(2.16) sup ||w|| L 2<T- sup ||6(-,t)|| L 2exp I - / Wd^u^Woo dr J . 

te[o,T] te[o,T] V^ Jo / 

In order to find a time interval such that the contraction property ( 12. 13ft holds, we need 
to estimate ||9 a: a(Mi)|| 00 and ||6(-,£)|| L 2. 
First we have 

lltZjjG^'UlJ Hoo W-Tuu'U'lx ~r " uu^lx\\oo 

< k{\\Ulx\\oo+ ||M1x||oo) 

< kfaR + aRWKWLi) 
=: C X R. 

The first term in (12.151) is bounded as 

(2.18) IIM^i) - a(u 2 )}v2 X \\ L 2 < CiR\\ui - u 2 \\ L ^ 

This can be seen from the following calculation: 
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[a(ui) - a{u 2 )}v 2x \\ L 2 = WiFuiui.Ux) - F u (u 2 ,U2)}v 2x \\ L 2 

< ||{F u (mi,Mi) - F u (u 2 ,Ui)}v 2x \\l^ + \\{F u (u 2 ,Ui) - F u (u 2 ,u 2 )}v 2x \\ L 2 



< CiRkl \\ui — u 2 \\ L 2 + \\ui — U2\\i 

< kciR(l + \\K\\ L i)\\ui - u 2 \\ L 2. 

If we assume F a (0, •) = 0, then the last term in ( 12. 15ft has a similar bound: 

(2.19) \\b{u 1 )-b{u2)\\ifl<C 2 R\\u\\r>. 
To obtain this bound, we decompose it the following way 

b(ui)-b(u 2 ) = -F a (u 1 ,u 1 ){u lx -U2 X }-u 2x {Fu(u 1 ,Ui)-F u (u2, ui)}-u 2x {F n {u 2 , Ui)-F c (u 2 , u 2 )} 
If we assume F a (0, •) = 0, we have F a (ui,Ui) = F au (£,,ui)ui, 

\\F a (ui,Ui){Uix -u 2x }\\ L 2 < k\\ui{ui x - U 2x }\\ L 2 

< kc R\\u lx -u 2x \\ L 2 

< kc R\\K x \\ L i\\ui -u 2 \\ L 2. 

Applying the mean value property to the remaining terms gives that 

||6(ui) -b{u 2 )\\ L 2 < k{c \\K x \\ L i +d\\K\\ L i +c 1 ||A'||^ 1 }i?||ni -u 2 || L 2. 
Substituting (lTTTD . (LTLSD and (I2TT9D into (I2TT6D . we obtain 

(2.20) sup \\v\\ L 2 < (d + C 2 )R-Te ClR - T sup ||w 1 -« 2 ||£3, 

te[o,T] te[o,T] 

which ensures (I2.13P if T < T 2 with 

T 2 



2e -(Ci + C 2 )R " 
Therefore, for < T < T* with 

T* = mm{T h T 2 } = — ^ min{l,i?}, 

the map S is a contraction on B^ in L°°([0,T];L x ) norm and thus possesses a unique 
fixed point u which is the unique solution of (II. ip . 

Note that without assuming F a (0, •) = 0, a different bound than (I2.19P is obtained 

||%i)-%2)||l 2 <(C4 + C 3 A)||£|| l2 , 
hence T 2 satisfying 

T 2 < 



2e{(C 1 + C 2 )R + C 4 } 

still ensures the contraction. This ends the existence proof. 

We prove the second part of Theorem 11.11 through the following corollary: 



u x \\ooYdr), 0<t<T 
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Corollary 2.4. Let u be the solution obtained in Theorem \l.l\ with a maximum life span 
[0,T). Then 

(2.21) \\u{t r )\\ H 2 < KH^exp (k(l + Cl )(l + \\K\\ w i,i) 4 f (1- 

where c\ is the embedding constant. This infers that only one of the following occurs 
i) T = oo and u is a global solution; 
ii) < T < oo and 

lim \\d x u(t, -)||l°° = oo. 

Proof. We use again the estimate in (I2.1(jp . setting v = u, 

II II ^ W II II II II II II II II u\\ 

-r;||M||^2 < — 1 1 a^ 1 1 oo 1 1 tt 1 1 ja-2 + 1 1 War 1 1 oo 1 1 «cccc 1 1 i 2 + \\ b \\H 2 - 

From a x = F uu u x + F uQ u x , it follows that Ha^oo < C\k{l + \\K\\ L i)\\u\\ H 2. Together with 
the estimates of Ha^H^ and ||6||h 2 m (12.51) and (12.31) . respectively, we obtain 

^-\\u\\ H2 < fc(l + Cl )(l + |H|oo) 2 (l + ||^||l0 3 (1 + ll^.|| L i)||n|| H2 . 
at 

Upon integration, we obtain (12.211) . The claim in ii) follows from a contradiction argument: 

If \im t ^T- || Wa | |oo < oo? it would lead to the boundedness of ||m||h 2 - One may therefore 

extend the solution for some T > T, which contradicts the assumption that T < oo is a 

maximal existence interval. □ 

3. Sub-thresholds for finite time shock formation 

3.1. Proof of Theorem 11.21 In this subsection, we consider the traffic flow model with 
Arrhenius look-ahead dynamics: 

d t u + d x (u(l - u)e- u ) = 0, 



w(0, x) = u (x) 



where u(t,x) = - j x 1 u(t,y)dy. Here 7 > denotes look-ahead distance. In the theory 



7 Jx 



of traffic flow, u(t,x) represents a vehicle density normalized in the interval [0, 1]. 

We want identify some threshold condition for the shock formation of solutions to 
(13. ip . From Corollary 12.41 we know that it suffices to track the dynamics of u x . Our idea 
is based on tracing M(t) := sup xeR [u x (x, t)] and N(t) := mi x£ ^[u x (x, t)}. The existence 
and differentiability (in almost everywhere sense) of M(t) and N(t) are proved in pQ. 

We also state a useful result, which is proved in [TT] . 

Lemma 3.1. (Lemma 3.1. in [TTJ^) Consider the following quadratic equality for A(t) 

d A 
(3.2) * =amA -b 1 (t))(A-b 2 (t)), A(0) = A , 

at 

with a(t) > 0, bi(t) < b 2 {t) and that a(t), b\{t), 6 2 (t) are uniformly bounded, 
i) If Aq > max 62, then A(t) will experience a finite time blow-up. 
ii) If there exists a constant b such that 

h(t)<b<b 2 (t), 
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then ( 13. 2\\ admits a unique global bounded solution satisfying 

min {Aq, min bi} < A(t) < b, 
provided Aq <b. 

With this result we obtain the following: 
Lemma 3.2. Consider the following quadratic inequality, 

(3.3) ^ > a(t)(B - 6i(t))(5 - b 2 (t)), B(0) = B , 

wj& a(t) > 0, b\(t) < 62(2) an ^ ^oi o(t), 61 (t), 62 (^) are uniformly bounded. 

i) If Bo > maxl)2, t/ien -B(£) will experience a finite time blow-up. 

ii) min{i?o,min6 1 } < B(t), fort > as long as B(t) remains finite on the time interval 

[0,t]. 

Proof, i) Subtracting (j3.2p from (J3.3P gives 

j(B -A)> a(t)(B - A)(S + A-h- b 2 ). 
Integration leads to 

(3.4) (B - A)(t) > (B - Ao) exp ( f a{t)(B + A-b 1 -b 2 ) dr 

Therefore, B > A implies B(t) > A(t). For any B > maxfe 2 set A = B , then by 
Lemma I3.1[ Aq will lead to a finite time blow-up of A(t). Hence, by (j3.4p . B(t) will 
experience a finite time blow-up. 

ii) Consider f l3.2|) . it is easy to see that min{Ao,min&i} < A(t). Then ( 13. 4p gives the 
result. □ 

We remark that Lemma 13.21 remains valid even if the quadratic inequality holds almost 
everywhere. 

Now, we are ready to prove Theorem 11.21 

Proof of Theorem 11.21 Let d := u x and apply dt to the first equation of (13. ip . 
d ■= (d t + (1 - 2u)e- a d x )d 

(3-5) _r , 

= e~ u 2d 2 + 2(1 - 2u)u x d - u(\ - u){u x } 2 + u(l - u)u x 

Define for t e [0,T), 

M(t) :=sup[u x (t,x)} =d(t,£(t)), 

(3.6) xm 

N(t):=w£[u x (t,x)] = d(t,7)(t))- 

The existence of £(£) and r)(t) is justified by Theorem 2.1 in [TJ. Furthermore, M(t) > 
and N(t) < 0. Then, along (£,£(£)), we have 

m„ = -K(e + 7) - «x(0} > -(~ M + Ar ). 

7 7 



12 YONGKI LEE AND HAILIANG LIU 

and (13.51) can be written as, 

M = e" s [2M 2 + 2(1 - 2u)u x M - w(l - u){u x } 2 + u(l - u)u xx 
(3.7) 

> e _fi ( 2M 2 + 2(1 - 2u)u ir M - u(l - u){u x .} 2 + u(l - it) ^~ M + N > 

And along (t,rj(t)), we have 

M ra = -{«x(»7 + 7) - U x (t])} > 0, 

7 
and (13.51) can be written as, 

N = e _a (2N 2 + 2(1 - 2u)u x N - u(l - w)^} 2 + u(l - u)u a 



(3.8) 

> e _0 ( 2iV 2 + 2(1 - 2u)u x N - u(l - m){m x } 2 



(I3.8P can be written as 

(3.9) N>2e~ ii (N-N 1 )(N-N 2 ), 
where 

iVl(M,Mx) 

and 

N 2 (u,u x ) 

We note that M < < N 2 because < u(t) < 1. It can be shown later that M is 
uniformly bounded from below, 

(3.10) JVi > --. 

7 

Applying Lemma [3.21 (ii) to ( 13. 9 I) with min 0<u<1 u<i Ni(u,oj) = — -, we obtain 



-(1- 


- 2u)u x - y/{(l - 2u)u x } 2 + 2«(1 - 


-«K 




2 




-(1- 


- 2w)w a . + y/{(l - 2u)u x } 2 + 2w(l - 


- M )^x 



7 



Mt) > mini - -, MO) 1 =: — . 

I 7 J 7 

Substituting this lower bound into ( 13. 7J , we obtain 

M > e- fi (2M 2 + J2(l - 2u)u x - U{1 ~ U) \ M - u(l - n)n 2 + " (1 ~ 2 ^° 

Rewriting of this inequality gives 

(3.11) M >2e~ il (M-M 1 )(M -M 2 ), 

where M 2 (> Mi) is given by 



"(!-") 1 1 . //on _ 9,,^ _ »(!-») 12 I 0„,n _ ,,^"72 _ Q^(l-n)ATo 



-{2(1 - 2u)u x - ^} + ^{2(1 - 2u)u x - U -^Y + 8n(l - tQg - 
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We claim that M 2 has an uniform upper bound, 



(3.12) 

By Lemma [3.21 (i), if 



M 2 < 



7 



1 y/2 



3-iVo 



M (0) > - 

7 



5 + ^-Vs-JV. 



then M(t) will blow up a finite time. This is exactly the threshold condition as stated in 
Theorem O 

To complete our proof we still need to verify both claims (j3.12p and (13.10p . 



To verify (13.121) . we set 

v :— 7 • u x — u(x + 7) — u(x). 
From < u(t) < 1 it follows that — 1 < v < 1. If suffices to find upper bound for M 2 over 



the set 
In fact, 



tt := {(u,v) e R 2 I < u < 1, -1 < v < 1}. 



{2(1 - 2u)t> - u(l - u)} + Vl 2 ! 1 - 2u ) v ~ «(1 - M )} 2 + 8u(l - ^)(« 2 - No) 

47 



<_[2+^4 + 2(l-JV )]. 

Here, we use maxr Uj v)e£i{~ 2(1 — 2u)v + u(l — u)} = 2 which can be verified easily since the 
underlying function is linear in v and quadratic in u. For the next one, max( n „) e n{8n(l — 
u)(v 2 — Nq)} = 2(1 —Nq) is used, which is obtained from the upper bound u(l — u) < 1/4. 
Finally, we are left with the verification of (I3.10p . With v defined above, we have 



Q := 1 N 1 



■(1 - 2u)v - y/{(l-2u)v} 2 + 2u(l-u)i 



By rearranging, 
(3.13) 

It follows that 

(3.14) 



a 1 



u(l — u)v 2 



Q ■ (1 - 2u)v 



< 



u(l — u)v z 



eQ s 



[l-2u)\. 



4e 



v\ 0<e<l. 



;i - c)Q 2 < ^{(1 - 2uf + 2ew(l - u)} 



< 



4e' 
where the maximum value is achieved at dVL. This gives 

Q 2 < 



46(1 - 6) 



Since e is arbitrary, we choose e = | to get Q 2 < 1, hence Q > — 1, which gives (I3.10p . 
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3.2. Proof of Theorem 11.31 We rewrite the traffic flow model ( 11. 3ft with the linear 
potential as 

(3.15) d t u + d x (u{l - u)e~^) = 0, 
where 

(3.16) u(t,x) = - (l + ^^)u(t,y)dy. 

IJx 7 

Let d := u x and apply c^ to ( 13. 15ft . 

d = (fifc + (1 - 2u)e^d x )d 

2d 2 + 2(1 - 2u)u x ,<i - u(l - u)-^.} 2 + u(l - u)u a 



(3.17) 
Here, 



2 f 1 /' x+7 
u x = <^ u(x) / u(y) dy 

7 I 7 Jx 

7 



7 



■(«-«), 



(3.18) 

where u— - J x 7 w(y) A/ as defined in the previous section. Define for t G [0, T), 



(3.19) 



M(t) :=sup{u x (t,x)] = d(t,£(t)), 

x£R 

N(t) := mi[u x (t,x)] = d(t,T)(t)). 



The existence of £(£) and ?](£) is justified by Theorem 2.1 in pp. Then, along (£,£(£)), 
( I3.17P can be written as, 



(3.20) 



M = e-*[ 2M 2 + 2(1 - 2u)u x M - u{\ - u){u x } 2 + u{\ - u)u x 



> e" a 2M 2 + 2(1 - 2u)u x M - u{\ - u){u x } 2 + u{\ - u) 



2{N - M) 

7 



where the last inequality follows from the fact that 

u xx (t, = -(u x - M) > -(N - M). 

7 7 

And along (t,rj(t)), ( 13.17ft can be written as, 

A> = e~" ( 2N 2 + 2(1 - 2u)u x N - u{\ - u){u x } 2 + w(l - u)u 2 

(3.21) ) 
> e" fi ( 2iV 2 + 2(1 - 2w)w x iV - u(l - w){w x } 2 

where the last inequality follows from the fact that u xx (t, rj) = 2 (u x — N) > 0. (13.2ip can 
be written as 

(3.22) N>2e~ il (N-N 1 )(N-N2), 
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where 



and 



-(1 - 2u)u x - y/{(l - 2u)u x } 2 + 2«(1 - 


-«)«x 


2 


-(1 - 2u)u x + y/{{\ - 2u)u x } 2 + 2m(1 - 


-M)u2 



iVi 



No 



We note that JVi < < N 2 because < u(t) < 1. 

By using the fact that < u < 1, and —2 < 7^3. < 2, it can be shown that Ni is 
uniformly bounded from below, 

iVi > — . 

7 



The verification of this inequality is similar to the one in the proof (I3.10p . details are 
omitted. With the lower bound of N\(t), Lemma F3.2I (ii) when applied to (13.221) gives 



(3.23) 



N(t) > min <^ - -, N(0) 

7 



iVo 

7 ' 



Substituting this lower bound into f l3.20|) . we obtain 

M>e" s 

= 2e~ u (M - Mi)(M - M 2 ). 



(3.24) 



2u(l-u)) , „ 2u(l -u)N 

2M 2 + \ 2(1 - 2m)m x i- -^ !>M - u(l - u)ul ' 



7 



T 



In order to apply Lemma [3721 (i) to (I3.24|) . we proceed to find the upper bound of M 2 (> 
Mi). Let v :— 7 • u x = — 2(m — w), then from the fact that < u, m < 1, we know that 
-2 < v < 2. We also let 

Q ■= {( u , v) G R 2 I < u < 1, -2 < u < 2} 

then M2 and it's upper bound are given by 

(3.25) 



M, 



{2(1 - 2u)v - 2u(l - u)} + J {2(1 - 2n)w - 2n(l - n)} 2 + 8u(l - u)(v 2 - 2N ) 



h 



< 



47 



16 + 2(4 - 2N 



Here, we use max( v ) e! ]{- 2(1 — 2u)v + w(l — m)} = 4 which can be verified easily since 
the underlying function is linear in v and quadratic in u. We also use u(l — u) < | in 
bounding the term 8w(l — w)(f 2 — 2iVo). Therefore, by Lemma F3. 2 1 (i), if 



M(0) > - 

7 



l + --V6-iVo 



then M(t) experience a finite time blow up. Hence we obtain the desired result. 
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3.3. Proof of Theorem 11.41 We only sketch the proof since it is entirely similar to that 
in the previous sections. Let d := u x and apply d x to the first equation of (11.11) to obtain 

(3.26) (d t + F u ■ d x )d = -F uu d 2 - 2F uU u x d - F aa u 2 x - F u u xx . 

It can be shown that < u < m, and therefore 

\u\ < m\\K\\wi,i, \u x \ < m\\K\\wi,i. 

To find the bound of u xx , we define for t G [0, T), 

M(t) :=snp[u x (t,x)] = d(t^(t)), 

N(t):=\ni[u x {t,x)]=d{t,T 1 (t)). 

From ([L2]), it follows that 

u xx (t,x)= / K'(z)u x (t,x - z)dz - K(0)u x (t,x). 

J — oo 

Therefore, along £(£), 

K{0)(N-M)<u xx <0, 
and (I3.26P is reduced to 

(3.28) M > -F UU M 2 - 2F ua u x M - F UQ u 2 x - F u K(0)(N - M). 

Also, along r](t), 

0<u xx <K{0){M-N). 

and (I3.26P is reduced to 

(3.29) N > -F UU N 2 - 2F uU u x N - F m u 2 x = -F UU (N - N^N - N 2 ), 
where 

Ar / _ \ r U yU x W [-T u fi i 'uu-t ' uu) U x 
N^U.Ux) = — . 

"uti 

From ( 13. 29ft we infer the lower bound of N(t) as 

N(t) > min{iV(0), min Nl(u,v)} -. N , 

0<M<m,|ti|<m|[A"|| 1/l/ i i i 

Substituting this lower bound into (13.281) . we obtain 

M > -F UU M 2 - 2F uU u x M - F au u 2 x - F U K{0){N - M) 
= -F uu (M-M 1 )(M-M 2 ), 
where 



2F ua u x - F u K(0) + \/{2F uU u x - F a K(0)} 2 - A{F uu F m v? x + F uu F u K(0)N } 
M 2 (u,u x ) = 



-2F 



II II 



Therefore, by Lemma [3.21 (i), if 

M(0) > max M 2 {u,v) =: A(iV(0)), 

0<u<m,\v\<m\\K \\ w i ,1 

then M(t) will blow up in finite time. Hence we obtain the desired result. 
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